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Abstract. Clouds seem like an every-day experience. But - do we know how clouds form on brown 
dwarfs and extra-solar planets? How do they look like? Can we see them? What are they composed 
pi |l of? Cloud formation is an old-fashioned but still outstanding problem for the Earth atmosphere, 

D ■ and it has turned into a challenge for the modelling of brown dwarf and exo-planetary atmospheres. 

' Cloud formation imposes strong feedbacks on the atmospheric structure, not only due to the clouds 

\^ , own opacity, but also due to the depletion of the gas phase, possibly leaving behind a dynamic 

■ and still supersaturated atmosphere. I summarise the different approaches taken to model cloud 

formation in substellar atmospheres and workout their differences. Focusing on the phase-non- 
equilibrium approach to cloud formation, I demonstrate the inside we gain from detailed micro- 
physical modelling on for instance the material composition and grain size distribution inside the 
cloud layer on a Brown Dwarf atmosphere. A comparison study on four different cloud approaches 
in Brown Dwarf atmosphere simulations demonstrates possible uncertainties in interpretation of 
' observational data. 
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^ INTRODUCTION 

: 

Brown Dwarfs are object much more massive than planets, but they can become as 
cold as planets during their lifetime. Young gas giants and also so-called close-in giant 
gas planets can be as hot as Brown Dwarfs (e.g. iflfjll ). Brown Dwarfs can only 
nucleosynthesis energy during a short period of their lift time, which is enough to bust 
their luminosity to a level where they can be much easier observed than any planet 
at a similar distance. The largest similarity between Brown Dwarfs and planets is the 
rich, multi-phase atmospheric chemistry. Descending the main sequence makes Brown 
Dwarfs and possibly late-type M-dwarfs the only objects, that form like stars and that 
contain dust (i.e. small solid particles) in their atmosphere. After this discovery Q23Q , it 
became immediately clear that Brown Dwarfs were not the perfect example for a static 
atmosphere but feedback of the dust clouds on the atmospheric structure, the chemistry, 
and the radiative transfer would need attention. Other places of efficient dust formation 
are the circulstellar envelopes of AGB stars which are extremely dynamic systems due to 
the presence of dust (see review by A.C. Andersen, and contributions by e.g. S.Hofner, 
N. Thureau this volume). [IjJ have presented spectra for L-type Brown Dwarfs showing 
broad absorption features between 8/im and 1 1/im which is the wavelength range were 
silicate dust would absorb. Although the idea of dust absorbing in the atmosphere is 
thermodynamically correct, it is hard to proof its existence directly since we are basically 
searching for the absorption signature of a thin transparent layer on top of an optically 
thick wall ( lfl9h ). [20] concluded that inhomogeneities in cloud decks and the evolution 
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FIGURE 1. Interaction of multi-scale processes during cloud formation in Brown Dwarf atmo- 
spheres. Dust formation and turbulence act on small scales, and gravitational settling and element replen- 
ishment act on large scales. 



thereof can plausibly produce observed photometric /-band variations. Extensive studies 
searching for uncorrected, time-dependent variability due to a variable cloud coverage 
of the atmosphere, suggest a variability of 2-3% ([14]) or even 2-10% (Bill]) depending 
in the wavelength studied. Richardson et al. (2007) infer the presence of silicate haze 
based on secondary-eclipse Splitzer observations on giant gas-planets, a conclusion 
which was also put forward by rf29ll based on HST transmission spectra for HD 189733b 
(see also reviews by G. Tinetti and by J. Harrington, this volume). 

CLOUD FORMATION PROCESSES 

Three major complexes need to be modelled to account for cloud formation in model 
atmospheres and to produce observable informations. Figure \T\ demonstrates how con- 
densation, gravitational settling and element replenishment interact. 

• Condensation: 

Seed particle form from the gas phase via a net or a path of chemical reaction. Once 
the first solid surface appears from the gas phase, chemical surface reaction can 
proceed rapidly to grow a massive mantle such that the initially tiny seed growth to 
a macroscopic grain of /im-size. While the nucleation process is only possible in a 
highly supersaturated gas, the growth can take place if a solid is thermally stable. 
The counter-acting process is evaporation if the temperature becomes too high. The 
condensation process reduces the element abundances is the gas phase selectively 
regarding those element forming grain monomers (like MgSi03, T1O2, S1O2 in the 
MgSi0 3 [s]-, Ti0 2 [s]-, Si0 2 [s] -lattice etc.). 

• Gravitational settling 

Once the grains are present they start to fall into the atmosphere due to the strong 



gravitational attraction in Brown Dwarfs and in planets. Since the grain sizes are 
not constant but different grain sizes will appear with different grain numbers, the 
equation of motion would need to be solved for a grain size distribution. An easier 
case is to solve the equation of motion for a height-dependent mean grain size 
derived from a grain size distribution. 
Element replenishment 

A static atmosphere would be dust free ( ll4lll ). 1271. 1361 and |25] infer from their 



observations that convection mixes up (elementally rich) gas from inner and hotter 
atmospheric regions, a mechanism essential for allowing the dust cloud to persist 
inside the atmospheres. The friction of the gas fluid is not very large so that its 
inertia drives it forward once it has been set into motion. Such fluids are called 
turbulent where whirly motions exist on a large range of scale. Convection is just 
one of them, but the largest on which the turbulent fluid field is energetically 
generated. Hence, modelling element replenishment means to describe mixing, 
advection, overshooting, and turbulence. 

As soon as the first Brown Dwarf was discovered, model atmosphere predictions where 
needed, and three different ideas appeared regarding the phenomenology of cloud 
formation: 

i) Static: The cloud particle are so small that they remain where they formed. 
(Tsuji, Copper et al.), 

ii) Bottom-up: The cloud forms by up-mixing of uncondensed material. 
(Ackerman & Marley, Allard et al.) 

iii) Top-down: The cloud is formed by particles falling into the atmosphere and evolv- 
ing on their way downward. (Helling & Woitke) 

Table Q] summarises the cloud models and their variants presently applied in the litera- 
ture. The grain size is an opacity input quantity, and the very first models set it constant 
assuming the grains to hoover inside the atmosphere (Tsuji). Later, experiences from 
cloud observations on Earth led to a time-scale comparisons to derive the grain sizes 
(Allard et al.) or to an analytic parameterisation thereof (Cooper et al., Ackerman & 
Marley). Only very recently, a grain size distribution based on a model for dust forma- 
tion was derived (Helling & Woitke). The calculation of the dust cloud opacity needs 
further the material composition of the cloud. The modelling of this is very much linked 
to the model assumption regarding the formation process. Almost all models assume 
phase-equilibrium between the atmospheric gas and the cloud particles, hence the su- 
persaturation ratio S = 1, or just a 1% supers aturation. However, only the departure 
from this equilibrium state, i.e. S ^> 1, would allow the actual grain formation. As- 
suming phase-equilibrium allows to consider the thermal stability of a vast variety of 
potential condensates (e.g. IZfJ) which then can be considered to deplete the gas-phase. 
This phase-equilibrium approach bears the danger that one considers compounds which 
might never have a chance to actually form. 

Despite the variety of the assumptions made to model clouds in Brown Dwarf atmo- 
spheres (Tabled]), fundamental understanding has been added by considering the limit- 
ing cases case B / dusty and case C / cond by Tsuji and Allard et al., respectively. The 
dust is element sink and opacity source inside the atmospheres in case B / dusty which 



TABLE 1. Cloud models in Brown Dwarf atmosphere simulations. 
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represents the the mid-L-dwarf regime. Case C / cond do consider the dust as element 
sink only and neglect its opacity. This case represented the T-dwarfs regime very well. 
Hence, these two simple limiting cases have established our picture from L-dwarfs as 
being covered in thick clouds, and of the T-dwarfs as Brown Dwarfs with an elementally 
depleted atmosphere but without visible clouds. As question remains how the clouds 
disappear going from L- to T-dwarfs (L-T transition region, [24], [fTsh . The main idea 
is that the cloud thickness changes. This isparameterised by a critical temperature T CI 
([38]), or a sedimentation efficiency / se d (iflll. lfioll ). or the mixing efficiency (Allard et 
al.). As long as no consistent theory exists to allow the simulation of the L-T, such 
empirical parameterisations are needed and modellers try to constrain them from ob- 
servations ( 0711 : Stephens et al., this volume). This inconsistency complicates also the 
evolutionary modelling from the L- into the T-dwarf regime as discussed in Q34fl . 



THE PHASE-NON-EQUILIBRIUM APPROACH 

The cloud formation process can only proceed in non-equilibrium state, else the sys- 
tem would not change at all. A phase-transition will therefore require a phase-non- 
equilibrium, hence, the gas phase needs to be supersaturated (S ^> 1). Depending on 
time and grain size scales, the dust formation proceeds via 1 .) nucleation, followed by 2.) 
growth (Sect. ). Nucleation runs on much shorter time scales than the growth process and 
needs a highly supersaturated gas to allow a reaction efficiency high enough to counter- 
balance destructive backward reactions. The (not-yet) particles which subsequently form 
in the nucleation regime are macro-molecules and clusters (e.g. [I13L 12 81). The result is 
a nucleation rate which determines the number of dust particles. jlA \2l\\ have shown 
that TiC>2 is a suitable candidate for seed formation in oxygen-rich environments since 
enough molecular TiC>2 is available to allow an efficient homogeneous cluster forma- 
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FIGURE 2. Material composition (in volume fractions V s /Vtot> left) ar, d altitude-dependent grain size 
distributions ( f(a) in cm _3 cm _1 ; right) in a stationary dust cloud layer in a Brown Dwarf atmospheres. 



tion process, in fact much more efficient than for instance SiO- or Fe-nucleation. Once, 
such a first surface emerges from the gas phase, lots of other compounds are already 
thermally stable and can therefore grow by chemical surface reactions. Still, thermal 
stability is not the only criterion, but we need to know which chemical surface reactions 
are possible (see [18]). A straight forward observation from gas-phase chemical equi- 
librium calculations is, that only small molecules like MgO, MgOH, SiO, SiS, AI2O, 
Ca(OH)2 are abundant enough to allow high growth rates. This conclusion is supported 
by condensation experiments 

The treatment of Ti02-nucleation by a modified classical nucleation theory piQ . and 
the calculation of the grow process including the effect of drift in connection with a 
parametrised up-mixing of uncondensed gas ([42, IMHtli) allows to study the cloud's 
material composition, the altitude-dependent grain size distribution, and also the phase- 
equilibrium state of the dust-gas mixture in the cloud beside various other details of 
the dust complex like nucleation rate, growth velocity, dust-to-gas ratio, metallicity. 
Figure |2] (left) shows that the cloud deck is populated by silicate grains mainly made 
of MgSi03[s] and Mg2Si04[s] beside some iron-compounds. The MgSi03[s] drops 
with increasing temperate which increases the fraction of Si02[s] inside ps 1000K. All 
silicates are evaporated at T ~ 2000K, and the cloud particles are now composed of 
mainly Fe[s], and Al203[s] at even higher temperatures. Figure |2] (right) shows a number 
of selected grain size distribution function, f(a) [cm~ 3 cm -1 ], sampling the cloud layer 
of that model. The lowest curve (1) represent the outermost grain size distribution in the 
region of nucleation of those depicted here. All subsequent f(a) (2 - 7) move towards 
larger grain sizes (abscissa) but only numbers 1, 2, 3 show an increase in numbers of 
grains. 4, 5, 6 have achieved the maximum possible number of grains and only move in 
grain-space to the right, hence, they represent cloud layers with continuously increasing 
grain sizes. 4, 5, 6 also narrower which demonstrates that the cloud base is populates by 
a relatively narrow interval of grain sizes. These big grains are mainly made of Fe[s], or 
Al203[s] just before they evaporate. 

The phase-non-equilibrium approach discussed so far does not assume a certain 
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FIGURE 3. Effective supersaturation ratio (S e ff) of the dust compounds considered in Fig. [2] Note that 
S e ff has been defined in (Helling, Woitke & Thi 2008) to discuss the contribution of all growth reaction of 
a particular compound compared to the standard definition. 



supersaturation ratio, S, but calculates it according to the present thermodynamical 
conditions and element abundances. We are therefore in the position to study if and 
where the dust-gas mixture forming the cloud layer is in phase equilibrium or not. Hence, 
we can study where the cloud is thermally stable (5=1) and where is is effectively 
processed in a constructive (5 > 1) or destructive (5 < 1) way. Figure [3] depicts the 
effective supersaturation ratio, S e fiEL defined to take into account the contribution of all 
growth reaction for a particular compound. The compounds considered fall into two 
regimes: 

• rare-element compounds: 

Ti-, A1-, Ca-solids exhibit an extreme phase-non-equilibrium 

• abundant-element compounds: 

Mg-, Si-, Fe-solids achieve phase-equilibrium over some pressure scale height in 
the cloud once the seed particles serve as condensation surface 

Figure [3] shows that also Fe[s] remains in phase-non-equilibrium in almost the entire 
cloud layer, hence only the low-temperature condensates can achieve phase-equilibrium. 



1 The classical supersaturation ratio assumes that the compound monomer (like a Mg2Si04-molecule for 
Mg2SiO/|.[s]-solid) does exist in the gas phase. Since this is not always the case, the supersaturation ratio 
for every surface reaction needs to be calculated since here the gas-phase constituents are known. For 
further details, see lfl9ll . 




TABLE 2. Brown Dwarf model atmosphere codes in comparison. See also TableQ] 
s - solid, si - solid & liquid. 
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34 83 3 as opacity source s 

10 as element sinks s 

84 680 43 as opacity source si 

169 as element sinks si 

83 ~ 2200 5 as opacity source s 

~ 1700 as element sinks si 

40 338 7 as opacity source s 

7 as element sinks s 



Consequently, phase-non-equilibrium models would show larger gas-phase abundances 
regarding molecules containing element of thermally very stable compounds as demon- 
strated in lll7n . 



COMPARING DUST CLOUD MODELS 



"Where stability and equilibrium are inherent signs if an unchanging and 'static' world, 
the formation and evolution of complex physical systems is intimately coupled to the 
existence and stabilisation of non-equilibrium states of matter." [28] . 

Section already hints that present cloud models are well described by the above 
statement taken from Q28[|. Namely, one kind of approaches studies the possible end- 
state of the cloud formation processes which is the (phase-)equilibrium state, the other 
kind has its emphasis on the formation process which is the (phase-)non-equilibrium 
state. Both aim to understand the influence of atmospheric dust clouds on observable 
quantities by means of model atmospheres simulations. A comparison study of such 
models with emphasis on the dust approaches was initiated at a workshop in Leiden 
2006(H). Table [2] demonstrates that not only the cloud models are different (Table Q]) 
but also the way the results of the cloud model approaches are used in the radiative 
equilibrium calculation differ. The number of compounds considered as opacity species 
is usually smaller than the number of compounds used as element sinks, and the element 
abundance data are taken from different sources. Part of the reason are missing refractory 
index data and numerical difficulties in handling such complex absorbing systems. The 
element abundances are particular interesting since they determine the composition of 
the gas phase (beside the local temperature and density) from which the dust forms. 

Figure |4] compares the results from four different model atmospheres codes (Ta- 
ble [2]) for the stellar parameter combination of a typical solar- metallicity L-dwarf . The 



2 http://www.lorentzcenter.nl/lc/web/2006/203/info.php3?wsid=203 
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FIGURE 4. Test cases for complete atmospheric models for logg = 5.0, solar element abundance 
with T e ff = 1800K. Note: Different colours stand for different stellar atmospheres codes. Four models 
are plotted for the Tsuji-case (brown): long-short-dashed: T cl - =1700K (extended cloud), short-dashed: 
T cr =1800K, long-dashed: T cr =1900K (thin cloud), dotted: no dust opacity considered. Different line 
styles in log (a) indicate different homogeneous dust species in the Marley, Ackerman & Lodders-models. 



temperature-pressure structures, (T,p) (1 st panel), differ the most in the low-pressure 
regime above the cloud layer and in the (T,p) -range where the cloud sits. The models 
which produce a thick cloud layer show the backwarming effect due to the strong cloud 
opacity at rs 2000K. The dust-to-gas mass ratio, pd/p g (3 rd panel), demonstrates where 
most of the dust mass is located, and coincides with the backwarming effect. However, 
the pd/ p g maxima differ widely amongst the models. Such a strong difference must be 
attributed to the amplification of small differences in a coupled system such as a stel- 
lar atmosphere, because Helling et al. (2008) demonstrate much smaller differences in 
the Pd/pg maxima if the dust models alone are tested. The models also differ in the 
predicted mean grain sizes, (a) (4 th panel), which reflects the differences in the cloud 
model approaches the strongest amongst the quantities in Fig.|4] 



CONCLUSIONS 



The modelling of Brown Dwarf atmospheres has turned into an unexpected challenge 
as these atmosphere contain dust which causes strong feedbacks as opacity source and 
element sink. Different approaches are applied to model dust in Brown Dwarf atmo- 
spheres, and obeying its nature as coupled system of nonlinear equation, atmosphere 
simulations applying different cloud models do produce different results although the 
over-all cloud- structure appears robust. A fundamental difference amongst the models 
is the consideration of phase-equilibrium. Kinetic models did show that in particular 
compounds made or rare elements like Ti, Al, and Ca never achieve phase equilibrium, 
i.e. thermal stability, in the Brown Dwarf atmosphere. Furthermore, the transition from 
the L- into the T-dwarf regime is not understood in enough detail that consistent atmo- 
sphere simulations are available. This is a sever challenge for evolutionary models. 

What else needs to be done? A still outstanding problem is the modelling of turbu- 
lent dust formation in large scale simulations which is also of significance for planetary 
atmospheres or protoplanetary disks. The challenge is here not to lose the small-scale 
interactions between the chemistry and the fluid field which has been shown to support 
the turbulent nature of the fluid field and to even cause the appearance of medium-scale 
cloud structures. A second challenge is the possible coupling of the dust-ionised atmo- 
sphere to the object's magnetic field. [35] and ll20ll argue against an ionisation potential 
in Brown Dwarf atmospheres. However, drifting grains in a turbulent environment might 
make a re-consideration necessary. 
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